  module physics_const_mod
  
    use kinds_mod
    use mesh_const_mod
  
    implicit none
  
    real(r8), parameter :: earth_radius  = 6371.220d3 ! m
    real(r8), parameter :: earth_omega   = 7.292d-5   ! s-1
    real(r8), parameter :: earth_gravity = 9.80616d0  ! m^2 * s-2
        
    ! thermal dynamic constants
    real(r_kind),parameter :: Rd        = 287.0583
    real(r_kind),parameter :: Rv        = 461.5233
    real(r_kind),parameter :: p0        = 100000. ! base pressure in Pa
    real(r_kind),parameter :: Md        = 0.0289644
    real(r_kind),parameter :: Mv        = 0.01801528
    real(r_kind),parameter :: eqv       = Md / Mv - 1
    real(r_kind),parameter :: rqv       = eqv + 1
    real(r_kind),parameter :: eps       = Mv / Md
    real(r_kind),parameter :: Cpd       = 1004.5
    real(r_kind),parameter :: Cpv       = 1846.
    real(r_kind),parameter :: Cvd       = 717.5
    real(r_kind),parameter :: Cvv       = 1385.
    real(r_kind),parameter :: T0        = 273.15
    real(r_kind),parameter :: alpha_d   = Rd / Cpd
    real(r_kind),parameter :: gamma     = Cpd / Cvd ! 1.4
    real(r_kind),parameter :: T_reduce  = 0.0065    ! vertical lapse rate of temperature
    real(r_kind),parameter :: ALPHA_ROT = 0.0       ! the flow orientation angle
    
  end module physics_const_mod
